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ABSTRACT 

We present detailed models of low and intermediate-mass asymptotic giant branch (AGB) stars with and 
without the '**F(Q;,/?)^'Ne reaction included in the nuclear network, where the rate for this reaction has been re- 
cently experimentally evaluated for the first time. The lower and recommended measured rates for this reaction 
produce negligible changes to the stellar yields, whereas the upper limit of the rate affects the production of '^F 
and ^^Ne. The stellar yields increase by ~ 50% to up to a factor of 4.5 for '^F, and by factors of ~ 2 to 9.6 for 
^'Ne. While the '^F(Q:,/?)^'Ne reaction competes with '^O production, the extra protons released are captured 
by '^O to facilitate the '^0(/?,a)'^N(Q;,7)''F chain. The higher abundances of '''F obtained using the upper 
limit of the rate helps to match the [F/O] ratios observed in AGB stars, but only for large C/O ratios. Extra- 
mixing processes are proposed to help to solve this problem. Some evidence that the '^F(Q;,/7)^'Ne rate might 
be closer to its upper limit is provided by the fact that the higher calculated ^'Ne/^^Ne ratios in the He intershell 
provide an explanation for the Ne isotopic composition of silicon-carbide grains from AGB stars. This needs 
to be confirmed by future experiments of the '^F(a,p)^'Ne reaction rate. The availability of accurate fluorine 
yields from AGB stars will be fundamental for interpreting observations of this element in carbon-enhanced 
metal-poor stars. 

Subject headings: nuclear reactions, nucleosynthesis, abundances. Stars: AGB and post-AGB stars. Stars: 
Carbon, Stars: Population II 



1. INTRODUCTION 

Interest in flie '^F(a,p)^'Ne reaction (Q value=1.741 MeV) 
came from early pre-supernova (SN) models that suggested 
that the reaction might be important in the helium and car- 
bon burning regions during the SN. After the shock wave 
increases the internal temperature and density, the timescale 
for destruction of '^F via the (a,p) reaction is comparable to 
that of its /3'^-decay lifetime (Arnett & Truran 1969; Truran 
et al. 1978; Giesen 1987), where the laboratory half-life of 
'^F is Ti/2 = 109 minutes. The early work by Arnett & Truran 
(1969) used unpublished theoretical estimates from Fowler; 
these rates did not appear in Fowler et al. (1975), Harris et al. 
(1983) nor Caughlan & Fowler (1988), and are only valid for 
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r > 800 X lO^K. Until 2006 the only rate for the i^F(a,p)2iNe 
reaction was the theoretical estimate available in the Brus- 
sels nuclear reaction-rate library (Aikawa et al. 2005). The 
first experiment aimed at determining the '^F(Q;,p)^'Ne rate 
over a large range of stellar temperatures was carried out by 
Lee et al. (2007, in preparation). This experimental evalua- 
tion, when considering its associated uncertainties, presented 
significant differences compared to the theoretical rate, espe- 
cially at the low temperatures relevant for He-shell burning in 
AGB stars (T w 300 x 10''K). In this paper we investigate the 
effect of such differences on the nucleosynthesis occurring in 
AGB models of various initial mass and composition. 

These are stars of mass less than ^ 8Mq located in the high- 
luminosity and low-temperature region of the Hertzsprung- 
Russell diagram. They have evolved through core H and He 
burning, and are now sustained against gravitational collapse 
by alternate H and He-shell burning (see Herwig 2005, for a 
recent review). AGB stars are the site of nucleosynthesis and 
mixing processes that lead to the production of carbon, nitro- 
gen, fluorine and heavy elements such as barium and lead. 
The strong stellar winds associated with these stars ensure 
that the freshly synthesized material is expelled into the in- 
terstellar medium, making AGB stars major factories for the 
production of the elements in the Universe (Busso et al. 1999). 

The theoretical estimate of the "^F(a,p)^'Ne rate was 
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not present in our previous works (Lugaro et al. 2004; 
Karakas et al. 2006), although we had included the species 
'^F because of its important role in the reaction chain 
'''N(a,7)i^F(/J+i/)i^O leading to the production of ^^O in the 
He shell. In this note we include this reaction in the net- 
work and study its effects in detail because preliminary results 
showed an enhanced production of '^F when employing the 
new upper limit of the ^**F(Q;,p)^'Ne rate. This is of interest 
because AGB models do not synthesize enough '^F to match 
the [F/O] abundances observed in AGB stars (Jorissen et al. 
1992; Forestini et al. 1992). This negative result remains even 
after examining most of the current error bars of the many re- 
actions involved in the complex chain of production of '''p in 
AGB stars, such as the '4C(a,7)''^0 and the ^'^F(a,pf^Ne re- 
actions (Lugaro et al. 2004). There are still uncertainties in the 
stellar models that could affect the match to the observations, 
in particular extra-mixing processes, as proposed by Lugaro 
et al. (2004). However, we will not be able to accurately pin 
down the effects of such uncertain stellar processes while our 
estimates of the abundance of '^F in AGB stars are still un- 
dermined by uncertainties in the reaction rates involved. 

The cosmic origin of fluorine is not yet completely under- 
stood. Type II SN explosions (Woosley & Weaver 1995) and 
stellar winds from Wolf Rayet stars (Meynet & Arnould 2000) 
both play a significant role in producing this fragile element 
alongside AGB stars (Renda et al. 2004). Observationally, 
AGB stars and their progeny (e.g. post-AGB stars, planetary 
nebulae) are the only confirmed site of fluorine production 
thus far (Jorissen et al. 1992; Werner et al. 2005; Zhang & Liu 
2005; Pandey 2006), with no clear indication for enhanced F 
abundances resulting from the j^-process in a region shaped 
by past SNe (Federman et al. 2005). Moreover, the recent ob- 
servations of a greatly enhanced F abundance ([F/Fe] = 2.90) 
in a Carbon-Enhanced Metal-Poor (CEMP) halo star polluted 
via mass transfer from a companion during its AGB phase 
(Schuler et al. 2007) represents further strong motivation to 
better understand the details of the fluorine production mech- 
anism in AGB stars. 

The "^F(Q!,p)^'Ne reaction could also affect the abundance 
of ^^Ne in the He-shell of AGB stars. There is a long-standing 
puzzle concerning the isotopic composition of Ne measured in 
stellar silicon carbide (SiC) grains extracted from meteorites, 
which formed in the extended envelopes of carbon-rich AGB 
stars. About 40% of these grains contain ^^Ne and/or '^He of 
nucleosynthetic origin (Heck et al. 2007). Being a noble gas, 
Ne is believed to be ionized and implanted in the SiC dust 
during the very last phases of AGB evolution (Lewis et al. 
1994; Verchovsky et al. 2004). Measurements performed on 
a large number of grains show that the observed Ne compo- 
sition can be explained by the mixing of He-shell matter into 
the envelope material of AGB stars (Lewis et al. 1990; Gallino 
et al. 1990; Lewis et al. 1994; Heck et al. 2007). While the 
■^"^Ne/^^Ne ratios are well reproduced in this scenario^, the 
^'Ne/^^Ne ratios are higher than predicted by AGB models. 
Lewis et al. (1994) attributed the higher than predicted abun- 
dance of ^^Ne to spallation reactions where the grains are 
bombarded by cosmic rays during their residence time in the 
interstellar medium. These authors hence related the excesses 
of ^'Ne with respect to the values predicted by AGB models 

The extreme enrichment of ^"Ne in these materials is historically known 
as the Ne-E(H) component in meteorites, whose presence was one of the keys 
leading to the discovery of stellar SiC grains in meteorites (Anders & Ziimer 
1993). 
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Fig. 1. — Reaction rate of '*F(a,p)^'Ne including the upper and lower 
limits. Also shown is the Brussels theoretical estimate of this rate. In the 
lower panel, the ratios of the current upper and lower limits with respect to 
the Brussels rate are shown. 



to the age of the grains. However, Ott & Begemann (2000) 
have shown experimentally that the majority of pre-solar SiC 
grains would have essentially lost all the ^'Ne produced dur- 
ing spallation by recoil. These authors suggest that the ob- 
served variations of the ^^Ne/^^Ne ratios in SiC grains are 
more likely due to the effect of nucleosynthesis in the He- 
burning shell of their parent AGB star, and this is also indi- 
cated by their correlation with nucleosynthesis effects in the 
Kr isotopic ratios. The identification of such nucleosynthe- 
sis effects, however, are to date missing. The ^^F(a,p)^^Ne 
reaction could play a role in this puzzle. 

For these reasons we aim to explore in detail the effect 
of the new experimental evaluation of the '^F(a:,p)^'Ne rate, 
briefly described in §2, on the production of fluorine and ^'Ne 
in detailed AGB models. Our methods and models are pre- 
sented in §3, results in §4 and §5, and we finish with a discus- 
sion and conclusions. 

2. THE "^F(a,f)2'NE REACTION RATE 

The measurement of the ^^F(a,p)^'Ne reaction cross sec- 
tion is made difficult by the short half-life of '^F. Owing to 
the problems associated with the production of a long-lived 
'^F target or a high intensity '**F beam, the first study of this 
reaction was based on the measurement of the time-reversed 
^'Ne(p,a)'^F reaction at the Dynamitron Tandem Laboratory 
Bochum (Giesen 1987). The cross section measurements at 
higher proton energies (Ep > 3 MeV) were based on the di- 
rect spectroscopy of the emitted a particles, while the lower 
energy range was investigated using the activation method by 
analyzing the decay of ^^F. The results were not published be- 
cause the low energy data were affected by the strong beam- 
induced background from the "*0(p,«)"*F reaction. With the 
development of an intense ^^F beam at the Cyclotron Re- 
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search Center at the Universite de Louvain la Neuve, a direct 
measurement of the "^F(Q:,p)^'Ne reaction became possible 
and the reaction cross section was determined at higher en- 
ergies (Ea > 1.4 MeV corresponding to Ep > 3.1 MeV, Lee 
et al. 2006). The beam intensity, however, was not sufficient 
to extend these first measurements to energies of relevance for 
stellar He burning. In a complementary experiment therefore, 
the low energy range (Ep < 2.3 MeV) of ^'Ne(p,a)'**F was 
re-investigated at the 4 MV Van de Graaff accelerator at the 
University of Notre Dame using the activation method (Lee 
2006, Lee et al. 2007, in preparation). 

The cross section of the '**0(/?,n)'**F background reaction 
was measured independently over the entire energy range. 
The corresponding '**F activity was normalized to the abun- 
dance of O impurities in the target, and subtracted from 
the ^'Ne(/7, a) induced '^F activity. Based on these data, a 
reaction rate for '*'F(Q!,p)^'Ne was determined for the stellar 
temperature range 0.2 < T (GK) < 1 .0 of relevance for AGB 
star nucleosynthesis. The lower limit of the cross-section 
measurement is mainly determined by the statistical uncer- 
tainty of the activation data, while the upper limit is based on 
the uncertainty associated with the '^O induced background. 
The resulting cross-section data were analyzed in terms of the 
R-matrix theory. The analysis, however, was hampered by 
the lack of detailed information about the specific parame- 
ters of the observed resonances. Fig. 1 shows the reaction 
rate as a function of temperature based on these recent mea- 
surements. The solid black line indicates the recommended 
rate while the long-dashed line and the dotted lines show the 
upper and lower linuts, respectively. These Umits of the re- 
action rate correlate with the experimental uncertainties in 
the cross-section data, as well as with the uncertainties from 
nuclear structure information. Shown for comparison is the 
predicted Hauser-Feshbach rate as a gray dot (Aikawa et al. 
2005). The present recommended rate is in good agreement 
with the Hauser-Feshbach prediction in the characteristic tem- 
perature range of AGB stars. Details of the experiment and 
the reaction rate analysis will be discussed in a forthcoming 
paper (Lee et al. 2007, in preparation). 

3. THE STELLAR MODELS 

The numerical method we use has been previously de- 
scribed in detail (Lugaro et al. 2004; Karakas et al. 2006). 
Here we summarize the main points relevant for this study. 
We computed the stellar structure first using the Mt Sttomlo 
Stellar Structure code (Lattanzio 1986), and then performed 
post-processing on that structure to obtain abundances for 77 
species, most of which are not included in the small stellar- 
structure network. This technique is valid for studying reac- 
tions not directly related to the main energy generation, as 
they can be assumed to have no impact on the stellar struc- 
ture This is certainly the case for studying the effect of the 
^^F(a,p)^^Ne reaction on AGB nucleosynthesis. On top of in- 
cluding neutron-capture reaction rates from Bao et al. (2000) 
for nuclei from Ne to S, the main change to the nuclear net- 
work for this study is the addition of the ^^F(a,p)^'Ne reac- 
tion rate into the 77 species network. 

The stellar-structure models used for this study are sum- 
marized in Table 1, and have been previously discussed in 
detail in Karakas & Lattanzio (2007, and references therein). 
Owing to the fact that we found the ^^F(a,p)^^Ne reaction 
to affect the abundance of '^F we have concentrated on mod- 
els that produce the most of it i.e. M ~ 3Mq (Lugaro et al. 
2004). We also show results from a lower mass (I.QMq) 



and two intermediate-mass (5Mq) AGB stars for compari- 
son. Both the 5Mq models experience proton-capture nucle- 
osynthesis at the base of the convective envelope (hot bottom 
burning, HBB). The 3Mq, Z = 0.012 model was computed 
with the revised solar abundances from Asplund et al. (2005), 
whereas the Z = 0.02 models were computed with Anders & 
Grevesse (1989) abundances. The lower metallicity models 
were computed using Anders & Grevesse (1989) scaled-solar 
abundances. We also present a model for a 2Mq Z =0.0001 
([Fe/H] ~ -2.3) star, which is relevant to the above-mentioned 
recent observation of highly-enhanced fluorine in a halo star 
of similar metallicity. 

A partial mixing zone (PMZ) is required to produce a '^C 
pocket in the He-intersheU during the interpulse period. It is in 
the '^C pocket that neutrons are released by the '^C(a,n)'^0 
reaction (Gallino et al. 1998); in this study we artificially in- 
clude a PMZ of constant mass at the deepest extent of each 
third dredge-up (TDU) mixing episode in exactly the same 
way as described by Lugaro et al. (2004). We include a pocket 
of O.OO2M0 for all lower mass cases, and we include a pocket 
of 1 X 1O"^M0 into the 5Mq, Z = 0.02 model. Note that these 
choices result in a '^C pocket between 10% to 15% of the 
mass of the He-intersheU region. 

In Table 1 we present the initial mass and metalUcity, Z, 
the C, N and O solar abundances used in the structure model 
where AG89 refers to Anders & Grevesse (1989) and A05 
to Asplund et al. (2005), the mass of the partial mixing 
zone (PMZ), the total number of thermal pulses (TPs) com- 
puted, the maximum temperature in the He-shell, TJ^'*", the 
maximum temperature at the base of the convective enve- 
lope, ^bce^' the total mass mixed into the envelope by TDU 
episodes, Massdred, and the final envelope mass M^nv All data 
are in solar units, except the temperatures, which are in mil- 
lions of kelvins. We present some information about the light 
elements including the surface C/O and ^^C/'^C number ratios 
at the last computed time step. 

4. RESULTS 

In Table 2 we show results from the stellar models that em- 
ployed the recommended rate of the ^^F{a,p)^^Nc reaction. 
For each mass and Z value, we show the C, N and O abun- 
dances used in the structure model (as for Table 1), the mass 
of the PMZ used in the computation, the yield (>) of '^F, 
the production factor (/) of '^F, and the multiplication fac- 
tor (X) needed to obtain the upper Umit '^F yield from the 
recommended-rate yield. All yields are in solar masses, the 
production factors / and the multiplication factors are dimen- 
sionless quantities. The same information is also presented 
for ^'Ne for each model. We compute stellar yields by inte- 
grating the surface abundances lost in the wind over the steUar 
lifetime, normalized to the initial abundance in the wind (see 
for e.g. Karakas & Lattanzio 2007). The production factors 
are defined according to / = logjQ(Xend/-^initiai), where Xend is 
the mass fraction at the tip of the AGB and ^initial is the initial 
mass fraction. The yields from the recommended calculations 
are essentially the same as the yields obtained from models 
that employed the lower limit, adopted the Brussels theoreti- 
cal rate, or did not include the '**F(a,/?) reaction at all. 

From inspection of Table 2 we can see that employing the 
new upper limit of the "^F(a,p)^'Ne reaction results in a 
significant increase in the production of ^^F and ^^Ne. The 
change in the yield increases with decreasing metallicity, at 
a given mass, with the largest change found in the 5Mq, 
Z = 0.004 model where the ''F yield increased by a factor 
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of 4.5. The largest change in the ^'Ne yield is a factor of 9.6 
for the 3Mq, Z = 0.008 model. While we find large increases 
in the F yield for both the intermediate-mass AGB models, 
the absolute yields are significantly smaller than those from 
the lower mass objects; this is because ^'F is destroyed by 
HBB. For example the 5Mq, Z = 0.02 model produced 3 
times less '^F than the 3Mq, Z = 0.02 case, whereas the 5Mq, 
Z = 0.004 model produced about 40 times less '^F than the 
2.5Mq, Z = 0.004 model. From Table 2 we note that the PMZ 
had little effect on the production of '^F and •^^Ne in the 5Mq, 
Z = 0.02 model. 

While increases in the ^^Ne yield as a consequence of us- 
ing the upper limit of the "^F(q;, p) rate are larger than for '^F, 
the overall amount of this isotope produced by AGB stars re- 
mains small. This is reflected in the production factors that 
are / < 0.3 for aU models but the 2.5Mq, Z = 0.004 and 2Mq, 
Z = 0.0001 models, where the production factors are 0.45 dex 
and 2.08 dex, respectively. The increase at very low metal- 
licity might be significant for chemical evolution studies of 
the Ne isotopes. Overall however, we conclude that the rare 
isotope ^'Ne is not significantly produced in AGB stars, even 
when using the upper limit of the ^^F(a,p) reaction in the cal- 
culations. Most of this isotope in the Galaxy originates from 
Type II SN (Woosley & Weaver 1995; Tinnmes et al. 1995), 
although it would stiU be an interesting exercise to include 
our AGB yields into a chemical evolution model. The impact 
of the upper limit on ^'Ne production is more important for 
stellar SiC grains, this is discussed further in §5. 

In this section we did not discuss the surprising result that 
the '^F(a,p)^^Ne reaction affects the production of '^F in 
AGB stars. It is not intuitive why this should be the case, 
so in the next section we outline the mechanism responsible 
for the production of the extra fluorine. 

4.1. The ^'^F production mechanism 

The enhanced abundance of '^F may be explained by con- 
sidering the '**0(/7,a)'^N(Q;,7)'^F reaction chain. Including 
the '^F(a,p)^'Ne reaction reduces the abundance of ^^O be- 
cause it competes with ^^O production via the ^^F(/?^zv)^^0 
decay. However, the extra amount of protons from (a,p) en- 
hances the "^0(/?,q;)''N reaction rate, even though ^^O pro- 
duction has been deprived from the decay. In other words, 
the sum NisQ+Np (where A', is the abundance by number of 
nucleus remains constant, however, the product NisoNp, on 
which the number of '^O+p reactions depends, is maximized 
when NisQ is equal to Np. 

We can analytically analyze the effect of the extra pro- 
tons on the ^^F production in the He-shell. We simpUfy the 
>^0(p,a)i^N(a,7)'^F reaction chain to the i^O(p,a)>^N re- 
action. In a He-rich region, all ^"*N is converted to ^^F via the 
(a, 7) reaction; this either decays to "^O via the /3'^-decay with 
a branching ratio of / or makes extra protons via the {a,p) 
reaction with a branching ratio of 1 - /. Then, the number 
density of '^O, Nisq, is written as fNu^, and for protons, Np, 
as A';,„ + (1 - f)N\4j^, where Np„ is the original number density 
of protons without the inclusion of the '**F(a,p)^^Ne reaction, 
and Nu^ is the ''^N from the H-burning ashes. Then, the reac- 
tion rate of "^0(/?, a)'^N can be written as 

NmoNp{av)(j,,a)=fNi4t^[Np„ + (l - f)N>4^]{<7v)(p^a) (1) 

= /M4nA^p„ ((7v)(p.„) X [ 1 + ( 1 - m^^/Npil) 

Since fNii^Np„ {av)(j,^a) is the rate of ^^0(p,a)'^N without in- 
cluding the ^^F(q!,p) reaction, the term [1 + (1 -/)M4n/A'p(|] 
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Fig. 2. — Ne isotopic ratios observed in meteoridc SiC grains and predicted 
in tlie intersliell of our 3, Z =0.02 and 1 .9Mq , Z =0.008 models. The plot is a 
reproduction of Fig. 8 of Lewis et al. (1994), where we have added the model 
predictions. For each model we plot the Ne isotopic ratios in the He intershell 
at the end of each TP occurring when C/0>1 is satisfied in the envelope of 
the star. The crossed full symbols represent models computed without the 
'*'F(Q.p)^'Ne reaction rate, which give a constant result. The full symbols 
represent models run using the upper limit of the '^F(a,p)^'Ne reaction rate. 
Dotted lines connect the "Normal Ne" component of solar composition to 
the "He-shell Ne" component corresponding to the final compositions of the 
intershell for the 3Mq, Z =0.02 model. 



may be thought of as an "'^F enhancement factor". The over- 
all ^^F production increases as long as Nu-f^/Np„ > 1, and this 
condition is well satisfied in the He-burning shell. During the 
network calculation a realistic Nuj^/Np,, ~ 10^"; this ratio is 
large enough to explain the enhanced fluorine production in 
the stellar models. 

As possible sources of uncertainty we can ignore the other 
'^'F +a channels, that is the (a,n) and the (a, 7). Accord- 
ing to the Brussels theoretical estimate (Aikawa et al. 2005) 
the (a, 7) reaction is approximately two orders of magnitude 
slower at 0.3GK than the ia,p), whereas the (a,n) is 40 or- 
ders of magnitude slower. 

5. 2iNE IN METEORITIC SIC GRAINS 

To address the puzzle of the ^'Ne/^^Ne ratio in stellar 
SiC grains we have analyzed the effect of using the new 
^^F(Q;,p)^'Ne reaction rate on the ^'Ne abundance in the He 
intershell of AGB stars. The results are shown in Fig. 2 and 
compared to the SiC data from Lewis et al. (1994). The plot 
is a reproduction of Fig. 8 of Lewis et al. (1994) where we 
have added our new model predictions. The SiC data are de- 
rived from measurements on samples of grains in bulk, i.e. 
collections of a large number (^millions) of grains. Differ- 
ent symbols represent measurements done on collections of 
grains sampling different sizes, from 0.01 to 5 /xm, as de- 
scribed in the figure. Note that, since measurements in bulk 
are performed on millions of grains, they can be only used to 
derive the average properties of the parent stars of the grains. 

Each data point in Fig. 2 is interpreted as having been pro- 
duced by a mixture between the material initially present in 
the envelope of the star and the material mixed from the He 
intershell into the envelope by TDU. These two "ingredients" 
are referred in the plot as the "Normal Ne" and the "He-sheU 
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Ne" components, respectively. The Normal component is 
taken to have solar composition. The SiC grains show a com- 
position dominated by a "He-shell" component extremely en- 
hanced in ^^Ne with respect to solar, as it is the composition 
of the He intershell of AGB stars. However, it is clear that the 
data points do not lie on the straight mixing line between the 
two components (the dotted lines in the plot), which means 
that the "He-shell" component must be variable if we want to 
account for all the different measurements. 

Model predictions presented in the plot are for the 3Mq 
Z =0.02 and l.QM© Z =0.008 models. These models are the 
best within our sample listed in Table 1 to represent the par- 
ent stars of SiC grains. This is because they reach carbon-rich 
conditions toward the end of their evolution (a necessary con- 
dition for the formation of SiC) and have masses (between 1 .5 
and 3Mq), and metallicities (close to solar) in the range of the 
best candidate SiC parent star models (see e.g. Lugaro et al. 
1999, 2003, for a thorough discussion). 

When we compute our models using the recommended, 
lower hmit, or Brussels theoretical evaluation of the 
^^F(a,;7)^'Ne reaction rate, the results are equivalent to the 
models computed without the inclusion of this reaction, and 
they are the same as those presented by Gallino et al. (1990). 
The ^^Ne/^^Ne ratio in the intershell is constant ~ 0.0004 
and the rightward shift to higher ^^Ne/^^Ne ratios observed 
in the grains cannot be reproduced. Note that in this case the 
abundances of ^''Ne and Ne are barely modified in the in- 
tershell, in particular ^'Ne is destroyed by factors 5 to 50 in 
the H-burning ashes and restored to its original Solar System 
value by neutron-capture reactions on ^"Ne during the TPs, 
with neutrons released by the ^^Ne(a,n)^^Mg reaction. One 
the other hand, models computed with the upper limit of the 
'^F(a, /?)^'Ne reaction rate show an increase in the ^'Ne abun- 
dance, and hence in the ^^Ne/^^Ne ratio in the intershell of up 
to a factor of 6, which is the number needed to reach up to 
the most extreme data point observed at ^'Ne/^^Ne=0.0033. 
The predicted intershell ^^Ne/^^Ne ratio increases with pulse 
number and with the stellar mass because the temperature in- 
creases and the '**F(o:,p)^'Ne reaction becomes more efficient. 
The last computed TPs reached 302 and 278 x lO^K for tiie 3 
Mq and the 1.9 M© models, respectively. 

Another possible way of producing a higher abundance of 
■^'Ne in the He intershell is by increasing the neutron-capture 
cross section of ^"Ne. The value we use is 0.199 mbarn at 
30 keV, which is reconmiended by Bao et al. (2000) and cor- 
responds to the experimental estimate of Winters & Mack- 
lin (1988). A much higher value of 1.5 mbarn at 30 keV 
was previously suggested by Almeida & Kaeppeler (1983), 
in which case the final ^'Ne/^^Ne ratio in the intershell of our 
3M© Z =0.02 model is equal to 0.002. However, the data of 
Almeida & Kaeppeler (1983) have recently been re-analyzed 
(M. Heil, personal communication) resulting in a much lower 
cross section of 0.303 mbarn at 30 keV. With this latest eval- 
uation the final ^'Ne/^^Ne ratio in the intershell of our 3Mq 
Z =0.02 model reaches only 0.00073. We also checked that 
possible changes in the neutron capture cross section of ^^Ne 
itself, and the current uncertainties of the ^^Ne(Q;,n)^^Mg re- 
action rate (Karakas et al. 2006) do not lead to significant vari- 
ations in the abundance of this isotope. These considerations 
lead us to conclude that the '^F(a,p)^'Ne reaction rate being 
close to its upper limit would be a pronusing explanation for 
the ^'Ne/^^Ne ratios in SiC grains. 

Finally, we note that increasing ^'Ne/^^Ne ratios are corre- 
lated with increasing ^^Kr/^^Kr ratios measured in SiC grains 
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Fig. 3. — Comparison of fluorine abundances observed by Jorissen et al. 
(1992) and model predictions for selected stellar models: 3Mq with Z=0.02, 
0.012, and 0.008; and 1.9M0 with Z=0.008. All models include a PMZ 
of 0.002Mq . Predictions are normalized in such way that the initial 
abundance corresponds to the average F abundance observed in K and M 
stars (see Jorissen et al. 1992). Crossed MS and S symbols denote stars with 
large N excesses. Each symbol on the prediction lines represents a TDU 
episode. Solid lines represent calculations performed using no "*F(Q:,p)2lNe 
reaction, which are equivalent to using the current lower hmit, recommended 
value and Brussels library rate. Dotted lines are calculations performed using 
the current upper limit of the rate. 



(see Fig.5 of Ott & Begemann 2000). This correlation can 
be qualitatively matched by considering that both the ■^^Ne 
and **''Kr abundances in the intershell increase with increas- 
ing temperature. This is because **''Kr is produced via the 
branching point at ^^Kr during the high-neutron density flux 
produced by the the ^^Ne(a,np^Mg reaction during TPs (see 
e.g. Abia et al. 2001). Quantitatively, however, our models 
can only match the lowest observed ^''Kr/**^Kr. It remains to 
be seen if this mis-match can be attributed to uncertainties in 
the nuclear properties of the ^^Kr branching point, or to inter- 
sheU temperatures higher than those of our models during the 
late AGB or the post- AGB phases. Further work is needed to 
address this point. 

6. DISCUSSION 

In Fig. 3 we show the evolution of the surface [^^F/'*0] ra- 
tio as function of the C/O ratio for four AGB models, com- 
pared to the observations of fluorine-enhanced stars from 
Jorissen et al. (1992). The models are selected to best rep- 
resent the features of the observed stars, that is, stars with 
masses in the range 1.5 to 3Mq (Wallerstein & Knapp 1998), 
and with metallicities around Z = 0.01. Similarly to the 
yields, the final surface ['''F/'^O] ratios from these models 
are roughly 50 to ~ 140% higher when calculations are done 
using the new upper limit of the '^F(a,p)^^Ne reaction. 

From Fig. 3 we see that using the new upper limit of the 
'**F(Q;,/?)^'Ne reaction can result in a match between the stel- 
lar models and the stars with the highest observed '^F abun- 
dances, but only for the very high C/O ratios of ~ 4-5, found 
in the 3Mq, Z = 0.008 model. In the lower mass models and 
in the 3Mq of solar metallicity, the new upper limit does not 
result in a match between the predicted and observed [F/O] 
abundances. Lugaro et al. (2004) suggested that extra-mixing 
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processes in AGB stars may help to solve this problem by con- 
verting C into N, hence decreasing the C/0 ratio for a given 
'^F abundance. Further indication of this possibility is the fact 
that for any given C/O ratio MS and S stars with the higher 
^'F abundance also have N excesses, and lower '^C/'^C ra- 
tios than predicted by standard models (Abia & Isern 1997). 
Detailed studies of the possible effects of extra-mixing phe- 
nomena are required, and will have to analyze the impact of 
using the higher '^F abundance obtained using the upper limit 
of the ^^F(a,/7)^^Ne reaction rate. 

There are many uncertainties that affect AGB stellar mod- 
els including the treatment of convection and mass loss (see 
Herwig 2005, for a detailed discussion). One modeling un- 
certainty that might affect the results is that most of the stellar 
models did not lose all of their convective envelopes when 
the evolution sequences ended, that is, they did not leave the 
AGB track, and could, in principle, experience extra TPs and 
TDU episodes. This possibility is discussed in Karakas et al. 
(2007), where it was estimated that one more TP may occur 
for e.g. the 3Mq, Z = 0.012 model. We do not repeat this 
exercise here owing to the uncertainty of the efficiency of the 
TDU at small envelope masses (see discussion in Karakas & 
Lattanzio 2007), but note that more TDU episodes would fur- 
ther enrich the '^F and '•^C abundances at the stellar surface. 

Another modeling uncertainty that will affect our results is 
the choice of mass-loss rate during the AGB. We used the 
Vassiliadis & Wood (1993) mass-loss prescription that was 
empirically derived from Mira-type variables and might over- 
estimate mass loss for semiregulars, thus terminating the TP- 
AGB phase too early and hampering the formation of C stars 
at low masses. At solar metalUcity we do not form carbon-rich 
stars with initial masses below 2.5Mq (Karakas et al. 2002), 
whereas typical C-star initial masses are 2Mq (Claussen 
et al. 1987), although this result is somewhat model depen- 
dent (Abia et al. 2001; Kahane et al. 2000). Regardless, this 
observational result is in contradiction to our models, and is 
caused partly by our choice of mass loss, and also because 
we do not find efficient enough (or any) TDU in the low-mass 
AGB models of » Z©. Certainly, a different choice of mass 
loss would have a significant effect on the stellar structure and 
on the resulting F and ^'Ne yields. We address this point in 
Karakas et al. (2006) for intermediate-mass AGB stars where 
the yields of ^^Mg and ^''Mg changed by more than an or- 
der of magnitude by using the Reimers mass-loss rate on the 
AGB instead of Vassiliadis & Wood (1993); we speculate that 
we would expect similar changes to the yields of lower mass 
stars but more work is needed to address this important point. 
One final comment is that the 3Mq, Z = 0.02 model becomes 
a C-rich star at a total (current) mass of ^ 2.3M0. Given the 
uncertainties in deriving total masses of C stars this is not en- 
tirely out of the range of expected C-star masses. 

In Fig. 4 we show the bolometric luminosities plotted 
against effective temperature for two stellar models that be- 
come C-rich near the lip of the TP-AGB phase. Fig. 4 can 
be compared to Fig. 6 in Busso et al. (2007), with lumi- 
nosities and temperatures from a selection of AGB stellar 
models computed with the FRANEC code (Straniero et al. 
2003), plotted against bolometric luminosities derived from 
observations of C-rich stars (see also Guandalini et al. 2006; 
Whitelock et al. 2006). In comparison to the FRANEC mod- 
els, AGB models computed with the Monash stellar structure 
code cover a similar range of Ten from 3,200 K to 2,500 K 
as most of the carbon stars, and cover the observed range of 
bolometric luminosities. Similar to the FRANEC models, we 




4500 4000 3500 3000 2500 2000 
Effective Temperature (K) 

Fig. 4. — Bolometric luminosity, M(,oi, versus effective temperature, T^ff, 
for the 3M0, Z = 0.02 (black dots), and the 1.9Mq, Z = 0.008 (gray crosses) 
models during the TP-AGB phase. The large variation in Mboi and Teff shown 
in this diagram is caused by the change in these observables during the AGB 
lifecycle (i.e. thermal pulse - dredge-up - interpulse). 



cannot match the T^s's of the coolest stars with temperatures 
~ 2000 K. However, we must be cautious about making con- 
clusions from this comparison because we are showing the en- 
tire AGB evolutionary sequence, not just the sequences when 
the model stars have C/O > 1 at the surface. Secondly, we do 
not include a realistic treatment of low-temperature molecu- 
lar opacities but instead we approximated the opacity from 
CN, CO, H2O and TiO using the formulations prescribed by 
Bessell et al. (1989), and corrected by Chiosi et al. (1993). 
These are fits to the molecular opacities of Alexander (1975) 
and Alexander et al. (1983), and while they do include a de- 
pendence on envelope composition do not treat properly treat 
C-rich compositions (see, for example, Marigo 2002). 

Both Marigo (2002) and Busso et al. (2007) have outlined 
the importance of using realistic low-temperature molecular 
opacities in detailed AGB models. Future work will study the 
effect of carbon-rich molecular opacities on the stellar struc- 
ture and nucleosynthesis. Marigo (2002) found that the in- 
clusion of C-rich opacities truncated the TP-AGB evolution 
fairly quickly (that is, in a couple of TPs) once the C/O ra- 
tio exceeded unity. This is because the C-rich molecules that 
form under such conditions caused the star to become larger 
and cooler, and this in turn increased the mass-loss rate. One 
of us (Karakas, Wood & Campbell, in preparation) is cur- 
rently studying the effect of such opacities on detailed AGB 
models and noticed similar trends, in that the evolution ends 
before the C/O ratio exceeds values much larger than 2. 
One consequence of this is that we would no longer predict 
the large C/O ratios found in the 3Mq, Z = 0.008 model (see 
Fig 3), and, subsequently, values of [F/O] greater than ~1. 

Another exciting future opportunity is represented by the 
comparison of our models of very low metallicity, e.g. 
Z =0.0001 (|Fe/H] ^ -2.3), to observations of fluorine in 
CEMP stars, which likely achieved their chemical peculiar- 
ities from an AGB companion. Schuler et al. (2007) ob- 
served [F/Fe] = 2.90 in one such star with a [Fe/H] = -2.5. 
The 2Mq Z =0.0001 model reached a fluorine production fac- 
tor of ~ 3.59, which translates into a huge [F/Fe] = 3.63. 
This may be enough to explain the observations of Schuler 
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et al. (2007), although we need to consider dilution due to bi- 
nary mass transfer The fact that fluorine production shows a 
strong dependence on the initial stellar mass (see, for exam- 
ple. Fig. 1 of Lugaro et al. 2004), suggests that we may use 
the detailed model predictions along with observations of the 
lowest metallicity stars, to provide constraints on the proper- 
ties of the initial mass function in the early Universe (see, for 
example, Tumlinson 2007). 

7. CONCLUSIONS 

In conclusions, the comparison of our results to observa- 
tions of [F/O] in AGB stars, and to the Ne composition of SiC 
grains suggests that the values of the '^F(a,/7)^'Ne reaction 
rate may he closer to the current upper limit. More exper- 
imental data for this reaction at temperatures below 0.4 GK 
are, however, required to help verify this result. The result for 
F in AGB stars is less compelling than the results for Ne in 
SiC grains, owing to the fact that we cannot match the whole 
observed [F/O] range. To add to this problem is the need for 
some extra-mixing process to alter the C and N abundances 
while not destroying '"^F. Also, AGB modeling uncertainties 
(e.g. mass loss and molecular opacities) could dramatically 
affect the predictions of F yields and surface abundances, ren- 
dering any conclusions uncertain. 

The modeling uncertainties related to extra-mixing, the 
TDU and mass loss do not affect, however, the intershell com- 
positions of our stellar models and thus do not apply to the 
discussion of the Ne composition of stellar SiC grains. From 
Fig. 2 and the related discussion, we see that the measured Ne 
isotopic compositions could be explained by the upper limit 
for the '**F(a,;5)^'Ne reaction. This tantalizing result is also a 
more reliable hint that the reaction is indeed closer to its upper 
Umit than the comparison to F in AGB stars. However, fur- 
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ther work is required to test this scenario, including a detailed 
investigation into Kr nucleosynthesis in AGB stars. 

Finally, the larger stellar yields of '^F obtained using the 
upper limit of the ^^F(a,p)^^Ne reaction should be tested in 
a galactic chemical evolution model of the type presented by 
Renda et al. (2004). An AGB contribution to the production 
of ^'Ne may also be considered, given that the upper limit 
of the '^F(a,/?) rate results in a larger production of this rare 
Ne isotope. The observations of low F abundances in stars 
in the globular cluster u Centauri by Cunha et al. (2003), 
where other observations clearly indicate pollution by AGB 
stars (e.g. Stanford et al. 2007), are puzzling. Clearly further 
work is required to address the nucleosynthetic origin of this 
most interesting and fragile element. 
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TABLE 1 

Data and results from the stellar models, see the text in §3 for details. 



Mass 


Z 


CNO^' 


PMZ 


TPs 






Mass^ii-^.^i 




c/o 


L/ C 


3.0 


0.02 


AG89 


0.002 


26 


302 


6.75 


8.1 (-2) 


0.676 


1.40 


118 


5.0 


0.02 


AG89 





24 


352 


64.5 


5.0(-2) 


1.500 


0.77 


7.84 


5.0 


0.02 


AG89 


lE-4 
















3.0 


0.012 


A05 


0.002 


22 


307 


7.23 


9.2(-2) 


0.806 


2.47 


168 


1.9 


0.008 


AG89 


0.002 


17 


278 


3.29 


2.2(-2) 


0.222 


1.30 


138 


3.0 


0.008 


AG89 


0.002 


29 


319 


10.5 


2.1(-1) 


0.549 


5.00 


519 


2.5 


0.004 


AG89 


0.002 


28 


308 


7.33 


1.9(-1) 


0.685 


11.9 


1300 


5.0 


0.004 


AG89 





81 


377 


84.4 


2.2(-l) 


1.14! 


2.64 


11.0 


2.0 


O.OOOl 


,\G89 


().()()2 


26 


307 


9.00 


2.2(-lj 


0.040 


105 


2.25(+4) 



^Initial CNO abundances where "AG89" refers to Anders & Grevesse (1989) initial solar or scaled solar abundances, and "A05" refers to Asplund et al. (2005) 
solar abundances. 



TABLE 2 

Stellar yields of and ^^Ne from the AGB models. 



Mass 


Z 


CNO 


PMZ 


M^'Frec) 






K^'NCrec) 




X(2iNe) 


3.0 


0.02 


AG89 


0.002 


5.84(-6) 


0.684 


1.526 


1.25(-6) 


0.053 


4.423 


5.0 


0.02 


AG89 





1.83(-6) 


0.223 


1.632 


2.09(-6) 


0.050 


3.463 


5.0 


0.02 


AG89 


lE-4 


1.87(-6) 


0.227 


1.625 


2.16(-6) 


0.052 


3.327 


3.0 


0.012 


A05 


0.002 


5.66(-6) 


0.676 


1.736 


1.39(-6) 


0.102 


5.330 


1.9 


0.008 


AG89 


0.002 


9.35(-7) 


0.583 


1.178 


1.60(-7) 


0.032 


2.340 


3.0 


0.008 


AG89 


0.002 


1.71(-5) 


1.466 


2.407 


4.52(-6) 


0.340 


9.609 


2.5 


0.004 


AG89 


0.002 


1.33(-5) 


1.752 


2.061 


2.81(-6) 


0.456 


8.364 


5.0 


0.004 


AG89 





1.45(-7) 


0.104 


4.582 


-2.58(-6) 


-0.627 


-1.965 


2.0 


0.0001 


AG89 


0.002 


1.67(-5) 


3.589 


1.975 


3.23(-6) 


2.080 


8.551 



